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Abstract. The time-dependent transport through single-molecule magnets coupled 
to magnetic or non-magnetic electrodes is studied in the framework of the generalized 
master equation method. We investigate the transient regime induced by the 
periodic switching of the source and drain contacts. If the electrodes have opposite 
magnetizations the quantum turnstile operation allows the stepwise writing of 
intermediate excited states. In turn, the transient currents provide a way to read these 
states. Within our approach we take into account both the uniaxial and transverse 
anisotropy. The latter may induce additional quantum tunneling processes which affect 
the efficiency of the proposed read-and-write scheme. An equally weighted mixture of 
molecular spin states can be prepared if one of the electrodes is ferromagnetic. 
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1. Introduction 

The single-molecule magnets (SMMs) are foreseen as building blocks of organic 
spintronic devices [1, 2]. Such systems generally behave as magnetic cores with a 
large localized spin and display slow relaxation of magnetization at low temperature 
mostly due to presence of the anisotropy-induced magnetic barrier [3]. Similar to 
quantum dot physics, two-terminal steady-state transport measurements performed 
on SMMs revealed charging effects such as Coulomb blockade, sequential tunneling 
or negative differential conductance [4, 5]. In the spin sector, Kondo related features 
were observed experimentally [6, 7, 8, 9] and investigated theoretically [10, 11]. The 
exchange interaction between the local molecular moment and the delocalized spins 
tunneling through the molecular orbital might be exploited to control the quantum 
state of the local moment, i.e. to ’write’ and ’read’ its quantum state [12], 

On the experimental side, various techniques [13] are currently used to attach the 
orbitals (ligands) surrounding the molecular magnetic core to the source and drain 
probes. Unlike standard transport setups used in quantum dot devices, molecular 
electronics requires more careful handling of the contact regions. The difficult task 
of isolating a single molecule between source and drain electrodes is nowadays realized 
by using more advanced methods as electromigration, mechanically controlled break 
junctions [14] or spin polarized STM [15]. Recently several groups pushed even further 
these techniques and reported controlled time-dependent transport measurements for 
such SMMs when the contacts were switched on and off by varying the substrate- 
STM tip spacing [16, 17] or by bending break junctions [18], and transient currents 
arising when a molecular tail couples to an STM tip have been recorded [19]. SMMs 
have also been integrated into carbon nanotube transistors to serve as detectors for 
the nanomechanical motion due to the strong spin-phonon coupling [20, 21]. In two 
cornerstone experiments Vincent et al. [22] and Thiele et al. [23] detected the nuclear 
spin of a single Tb^"*" ion embedded in a SMM together with the Rabi oscillations. 

These promising experiments motivated us to investigate the transient transport 
properties of SMMs, with special emphasis on the regime when the couplings to the 
source/drain electrodes are switched on and off periodically. The transport regime we 
are interested in is similar to the so called turnstile pumping setup which represents a 
long-standing [24] asset of pumping or pump-and-probe experiments with quantum dots 
[25]. Along this periodic pumping, the source and drain tunneling barriers open and 
close consecutively, such that a single electron is transmitted across the sample; details 
on the turnstile operation will be presented in Section 3. To our best knowledge the 
transient and the turnstile regimes have not been theoretically investigated so far in the 
context of transport across SMMs. 

On the theoretical side the magnetic interactions in SMMs are described by effective 
giant spin Hamiltonians [3], mostly because of the large value of the localized magnetic 
moment. Using this description, in Ref. [26] the authors investigated the role of 
relaxation on inelastic charge and spin transport across a SMM weakly coupled to 
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metallic gates. In Ref. [27] transport across a SMM coupled to two ferromagnetic leads 
was investigated and it was found that the spin current across the SMM can reverse 
the localized spin if the leads have opposite spin polarizations. Memristive [28] and 
thermoelectric [29] properties of SMMs were also investigated. All these studies convey 
a similar message: The current can induce magnetic switching of the localized magnetic 
moment if the applied bias voltage exceeds the gap between the ground and excited 
states. 

In the absence of transverse anisotropy the effect of quantum tunneling of 
magnetization (QTM) is negligible and the full magnetic switching requires the transient 
occupation of all excited (intermediate) states with magnetic quantum numbers in the 
range [—S', S']. When present, the QTM might leave its hngerprint on the transport 
properties at resonant values of an applied magnetic held [30]. 

A complementary approach to transport properties of SMMs relies on density 
functional theory (DFT) [31, 32]. In the DFT approach the molecular structure and 
the contact regions are carefully taken into account, while the many-body correlations 
within the SMM are accounted for within various approximations. A detailed ab initio 
Hubbard many-body model for molecular magnets has been recently implemented [33] 
and allows the calculation of magnetic interactions. 

In the present work we investigate transient transport and turnstile pumping across 
a SMM embedded between magnetic and non-magnetic electrodes. As we are interested 
in the time dependent evolution of the currents and the accumulation of the charge and 
spin on the SMM we rely our investigation on the generalized master equation (GME) 
technique [34, 35]. Let us stress that, to capture the turnstile regime, one has to go 
beyond the steady-state rate-equation approach. We hnd that by setting the SMM 
in the quantum turnstile (QT) conhguration with ferromagnetic leads one can address 
two new issues which are relevant for the use of molecular states as magnetic qubits: 
(i) the one-by-one all-electrical writing and reading of excited molecular states with 
spin S — 1,.., —S -|- 1 (S' being the molecular spin of the initial ground state) and (ii) 
the controlled preparation of statistical mixtures of such intermediate states. So far, 
the stepwise magnetic switching protocols for excited molecular states that we propose 
here have not been investigated. In fact, previous studies (see e.g. Refs. [27, 26]) were 
focused only on investigating the full magnetic switching. The second issue was partially 
addressed by Tejeda et al. [36] some time ago. Their proposal concerns the preparation 
of equal weight superposition of states (e.g. |'0) = l/\/2{\S) -|- [S' — 1))) using at least 
two molecular clusters embedded in micro-SQUIDs. 

The rest of the paper is organized as follows: In Section 2 we present the theoretical 
framework by introducing the model Hamiltonian and giving a summary of the GME 
method. Section 3 presents the main results of our work while in Section 4 we give the 
conclusions. 
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2. Theoretical framework 

2.1. Model Hamiltonian 

The setup that we consider here consists of a SMM coupled to two external electrodes 
(see the sketch in Fig. 1). We investigate time dependent transport in the sequential 
tunneling regime, in which the electrons tunnel one by one from the left (source) 
electrode to the unoccupied molecular orbitals of SMM and then escape to the 
right (drain) electrode. In the present work we neglect all other possible transport 
mechanisms, such as the cotunneling processes which are responsible for the Kondo 
effect [10]. 

The model Hamiltonian contains several terms describing the SMM itself {Hm), 
the left {Hl) and right {Hr) electrodes and the time-dependent tunneling part HR- 

Hit) = Hm + Hl + Hr + Hrit). (1) 

In general, SMMs are characterized by a large spin S > 1/2. They do also present 
transverse anisotropy and an easy axis of magnetization [3]. Assuming that the relevant 
contribution to the transport comes only from the lowest unoccupied molecular orbital 
(LUMO) the SMM can be modelled by an effective, minimal Hamiltonian [26, 27] of the 
form: 

HM = en + Uh^h^ - Js-S- DSl + E{Sl - S^) - gfiRBSl. (2) 

In its simplest form the LUMO orbital consists of a single spinful interacting level with 
energy e, on-site Coulomb energy U and occupation h = h^, and is coupled to 

the localized spin S through an exchange interaction with a coupling strength J. The 
fourth and the hfth terms in Eq. (2) describe the easy-axis and transverse anisotropy 
with the corresponding constants denoted by D and E. For certain molecules [38] 
E D but this is not always true, as E can grow up to Zl/3 in other situations [3]. 
In the presence of an external magnetic held B pointing in the ^-direction, a Zeeman 
term is supplemented in Eq. (2), with g and fiR being the gyromagnetic factor and the 
Bohr magneton, respectively. In view of further discussions we single out the transverse 
anisotropy term and write Hm in Eq. (2) as 

Hm = HI, + E{SI-SI). (3) 

The reason behind this separation is that, contrary to Hm, H^ has an extra abelian U (1) 
symmetry generated by the ^-component of the total spin §1 = s^. Consequently, 

the eigenstates of Hf[ can be organized according to the eigenvalues m of §1 [39]. On 
the other hand the transverse anisotropy term does not commute with §1 and needs to 
be treated separately. We shall discuss in more detail the eigenstates of H^ in section 

2 . 2 . 

The source and drain electrodes are modeled as spin-polarized one-dimensional 
discrete chains which in the momentum space representation are described by the 
Hamiltonians: 

Ha — ^ ^ / dqa £qaa '^h^a'^qaO-^ 

a Jo 


a = {L, R}. 


( 4 ) 
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Figure 1. (Color online) A sketch of the SMM coupled to source and drain electrodes 
via time-dependent tunneling barriers. The tunneling amplitudes are controlled by the 
switching functions xl and XR- The turnstile operation consists in turning periodically 
on and off the contacts. Note that the left contact opens first, for the charging sequence, 
while the right contact couples later for the discharge/depletion sequence. The chemical 
potentials of the leads are chosen such that pl > Pr\ the bias eV = pl — Pr- to is 
some initial time. The leads are spin polarized in the so called antiparallel (AP) 
configuration. 

Both leads present an energy dispersion law of the form Sq^ = 2rcosg + Ao-, with 
r the effective hopping between the nearest neighbor sites in the leads and A^. a rigid- 
band spin splitting that describes the polarization of the leads. In Eq. (4) creates 
an electron with momentum q and spin a in the lead a = {L, R}. 

The last term in Eq. (1) describes the hybridization of the SMM with the contacts 

HTit)= / dqaXa{t){V^alaq^^ + h.c), ( 5 ) 

a=L,R a ^0 

where aj. creates an electron with spin a on the LUMO orbital and is the hopping 
amplitude of a tunneling process between the LUMO and the majority {a = -|-) and 
minority (a = —) electron states in the lead a. The coupling of the SMM to the contacts 
in the case of collinear magnetic conhguration, and in the absence of the switching 
protocol {Xa(t) = 1), is described by T" = 27r|l/^p p„o-(0), where Pao-(O) is the spin 
density of states at the Fermi surface for electrons in lead a. 

In view of further investigations we allow tunable spin polarizations in the leads 
and define P“ := (T® — r")/(r" -|- T"). 

Note that the Hamiltonian iLr(t) contains two time-dependent dimensionless 
functions Xa{t) (a = L,R) which simulate the switching of the contacts between the 
molecule and the leads. As we are interested in the turnstile pumping it is enough to 
consider them simply as rectangular periodic pulses (see Fig. 1) 
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In this section we discuss the energy spectrum and the organization of the eigenstates of 
the SMM Hamiltonian Hm introduced in Eq. (2). We shall start by discussing first the 
spectrum of When J = 0, the LUMO orbital gets decoupled from the local spin 
and the Hamiltonian has three U(l) symmetries generated by the local charge Q 
accumulated on the LUMO and by the 2 ;-components of the LUMO and local spins, s^ 
and Sz. Consequently, {Q, Sz, *^ 2 } provide the quantum numbers according to which the 
multiplets of the Hamiltonian are classihed. Notice that this low symmetry classification 
is valid for hnite magnetic helds. When B = 0, the Hamiltonian has a much 
higher symmetry, i.e. Uq{1) x SUs{2) x SU^{2) in the charge and spin sectors, but this 
situation shall not be discussed here, as we always assume a finite magnetic held. In 
the case we consider here, the classihcation of the states is rather trivial and we can 
simply denote the eigenstates as follows: | 0 , 0 ,S' 2 ), |l,t, 5'^), |1,4,, S'^), and | 2 , 0 ,S' 2 ), 
with Sz = —S, —S +1, S. In the presence of Coulomb interaction the double occupied 
states |2, 0, Sz) have an energy of the order ~ U, which is the largest energy scale in the 
problem, and in view of the discussion that follows, shall not contribute to transport. 
Therefore, to simplify the notations, it is enough to relabel the states and keep track of 
the Sz and Sz quantum numbers. In this new notation we have |0, 0, Sz) |0, Sz) and 
|l,cr, S'^) ^ |(T, S'z). 

A hnite exchange coupling J breaks the three U(l) symmetries down to Uq{1) x 
Ugt (1) generated by LUMO charge Q and the ^-component of the total spin Si. Still, the 
Hamiltonian Hfj can be diagonalized exactly and the states constructed in an analytical 
fashion in terms of the states introduced previously for J = 0, by using the Clebsch- 
Gordan construction [26]. Now the new states \Q, m) shall be classihed by the molecular 
charge Q, and by the 2 ;-component of the total spin, m. 

For m e [—S + 1/2, S' — 1/2] the single-particle states {Q = 1) are given by: 

\l,m)^ = l,m + l/2) + C^^\ t,m- 1/2), (6) 

and their associated eigenvalues read as: 

e - gi^sBrn + J-D + ^ ± A£{m), (7) 

where AS{m) = [D{D — J)m^ + (J/4)^(2S' -|- 1)^]^/^. The coefhcients in Eq (6) 
are simply the Clebsch-Gordan coefhcients. The states |0,S'^) are not ahected by the 
exchange coupling and one has |0,m) = |0, Sz). The corresponding eigenvalue is simply 
i’o,m = —DSl — guBBSz. The remaining Q = 1 states are |1, —S' —1/2) and |1, S'-1-1/2). 
For a vanishing magnetic held, B = 0, the states associated to ±m are degenerate and 
one has 

SfJB = 0 ) = = 0 ), ( 8 ) 

^0,m{B = 0) = So_rn{B = 0). (9) 

So far we have discussed how to construct the states and to compute the energy 
spectrum for H^. In the rest of this paragraph we shall address the role of the transverse 
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anisotropy term. The transverse anisotropy term {S'i - {S‘^ + Sf) does not 

commute with S'* and induces transitions [3] between the states of with the selection 
rule |m — m'l =2. 

As the molecular charge is a good quantum number even in the presence of the 
transverse anisotropy, the eigenstates of the total molecular Hamiltonian Hm can be 
classified by the molecular charge Q only. We shall label them where Q = {0,1} 

(states with molecular charge Q = 2 are disregarded) while u is an internal label that 
indexes the states within a multiplet. In the presence of the transverse anisotropy E, 
the ‘empty’ molecular states (EMS) can be written as: 

\7^o,u) u = l,..,2S + l, (10) 

m 

with m running over all allowed values in the range [—S, S']. For half integer S' and a 
small magnetic held, the transverse anisotropy plays a minor role in the mixing of the 
states as the transition amplitudes between the empty molecular states {|0,m)} 

are negligible (see also the discussion following Fig. 2.) 

In contrast, the transverse anisotropy couples the degenerate, one-particle states 
{Q = 1) with opposite m’s. The strongest mixing is expected for the pairs |1,1)^ and 
|1,—1)^ as the off-diagonal matrix element ^(1, is linear in E. Higher 

order mixing effects become important as the ratio E/D increases and one can generally 
write: 

1 '=!,.. ,2(25+1). (11) 

m s 

The eigenvalues Eq^^, of Hm and the coefhcients in Eqs.(lO) and (11) can be found only 
by numerical diagonalization. More details on the spectral properties and state mixing 
will be given in Section 3. 

Finally, we write down the matrix elements of the tunneling Hamiltonian with 
respect to the eigenstates of and derive the selection rules for molecular transitions 
due to electronic back-and-forth processes namely (A^i- = 1, = — 1): 

^(l,m|4|0,m') = (12) 

This equation describes the tunneling of one electron with spin a on the SMM orbital, 
when the number of electrons in the molecule increases by one while the total magnetic 
quantum number m changes by ±1/2. 

2.3. Generalized master equation approach 

The GME approach which we use to investigate the time-dependent transport relies on 
the partitioning approach [37]. More precisely, transient currents develop in the source 
and drain electrodes as they are contacted to the molecule at some initial instant. 
The leads are viewed as non-interacting particle reservoirs with chemical potentials /i^ 
(cK = L,R), and at equilibrium described by the Hamiltonian (4). This setting is suitable 
for perturbative calculations with respect to the lead-molecule couplings T^, and allows 
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us to compute transient currents in the presence of time-dependent modulation of the 
contacts, as in the turnstile regime. 

The GME method essentially provides the SMM reduced density operator (RDO) p 
dehned as the partial trace over the leads’ degrees of freedom p{t) = Trj-^j{W(t)}. Here 
W{t) is the density operator of the whole structure which solves the Liouville - von 
Neumann equation ihW{t) = [H(t), W{t)], and the trace is over Jei - the Fock space of 
the non-interacting electronic reservoirs. In the sequential tunneling regime considered 
here the master equation takes a rather compact form (for a full derivation see Ref. 
134]): 

^ = - {P,|p|} , (13) 

0.[p 1 = [-ffT(t). f (14) 

Jo 

where we introduced the “free” evolution operator of the disconnected system Ut = 
^-i(HM+HL+Hii)t/h equilibrium distribution of the leads pei [40]. The dissipative 

operator TJf collects all sequential tunneling processes from the switching instant to = 0 
to the current time t. We solve numerically Eq. (13) with respect to the fully interacting 
states of Hm and obtain the populations associated to a given state \p>q,u) as 

PqA^) = (15) 

Once p{t) is known one can calculate average values of molecular observables by 
performing a trace over the Fock space JFm of the molecule. For instance, the total 
charge accumulated on the orbital involved in transport is given by 

Q{t) = eTr^^{p(t)iV}, (16) 

where the total electronic occupation N = + fi]^ and e is the electron charge. The 

continuity equation then becomes [35] 

■hit) - JRit) = ^^{^Q,uQ\pit)WQ,uQ), (17) 

Q tdQ 

where z/g is the set of states with charge Q . By inserting the tunneling Hamiltonian 
Ht into the double commutator Vt given in Eq.(14) one identihes Jl and Jr from 
the RHS of Eq.(17). It is straightforward to show that the ‘empty’ molecular states 
\<Po,v) do not contribute to the currents. Similarly, one can calculate the total spin 
{S\) = TTRj^{p(t)Sl} as well as the spin currents. In this work the relaxation of 
the excited molecular states via phonon emission is not considered. This is a good 
approximation as long as the timescale on which the quantum turnstile operates is much 
smaller than the relaxation time which is of order of 10“® s (see e.g Refs. [41, 42]). In 
fact previous studies [43] reported that the current-induced magnetic switching is stable 
against intrinsic spin-relaxation processes. 
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Figure 2. (Color online) (a) The energy levels for a SMM of spin S = 5/2 as function 
of the total magnetic quantum number m. The two-particle sector is not included as 
the corresponding states are outside the selected bias window. Other parameters are: 
e = 0.25meV, J = O.lmeV, U = ImeV, D = O.OdmeV, ggsB = 0 and E = 0. (b) The 
transverse anisotropy E induces mixing of degenerate or nearly degenerate states with 
TO = 1 and TO = — 1 (indicated by the double arrow). There is no significant quantum 
tunneling of magnetization between the states |0 ,=f 1/2) and |0,±3/2) (indicated by 
the dashed lines). The numbers in the circles are the total magnetic quantum numbers 
TO. (c) Tunneling processes connecting the lowest four states of the SMM via back- 
and-forth tunneling with the leads. These states are also involved in the first turnstile 
cycle - see the discussion in the text. 


3. Results and discussion 

3.1. The transport configuration and tunneling processes 

The numerical simulations were performed for molecules with spin 5 = 5/2 but our 
conclusions remain valid for larger half-integer values of S. For reasons that will become 
clear below in this work we restrict ourselves to SMMs with small transverse anisotropy, 
that is E D J (in fact we allow a maximum ratio E/D = 1 /25 a.t hxed easy-axis 
anisotropy constant D). We shall investigate two spin conhgurations for the leads. In the 
so called antiparallel (AP) conhguration the left lead carries only spin-down electrons 
and the right lead is spin-up polarized. For the second conhguration the left lead is 
non-magnetic (i.e. = 0) and the right lead remains ferromagnetic. We label this 

conhguration as normal-ferromagnetic (NF). 

The chemical potential of the leads fiL,R are set such that only the states and 
contribute to the tunneling processes, while double occupied states, \(p 2 ,u)-i have a 
much higher energies, E 2 ^u > Tl, and do not contribute to transport. 

Let us hrst discuss the energy spectrum and the relevant lead-molecule tunneling 
processes in the absence of transverse anisotropy. In the following discussion we shall 
use the basis {\Q,m)} of = Hm{E = 0). As we are interested in the pumping 
mechanism at hnite transverse anisotropy E, we shall switch later to the basis \<PQ^y) of 
the full Hm- In that situation, with certain modihcations, a similar turnstile scenario 
holds. Figure 2(a) shows the energy levels of for a given set of parameters, and in 
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the absence of the external magnetic field. Figure 2(b) schematically shows the charge 
Q = 1 (integer m’s) and Q = 0 (half integer m’s) branches of the spectrum; in view 
of further discussion the double arrow and the dotted lines mark some of the quantum 
tunneling of magnetization (QTM) processes induced by a nonvanishing E. The states 
connected by the double arrow are strongly mixed by E, while the ones connected by 
the dashed lines are only weakly coupled. As charge Q is conserved, it implies that 
direct transitions between Q = 0 and Q = 1 branches are forbidden by symmetry. This 
is only possible through processes involving states in the leads that do not conserve the 
charge. For example in Fig. 2(c), we show such processes (blue and red arrows). The 
same hgure also shows how the SMM evolves from an initial ‘empty’ molecular state 
|0,5/2) to the to the next EMS |0,3/2) via tunneling processes. 

We shall call the transitions m m — 1/2 ‘forward’ processes (they follow the 
full line arrows in Fig. 2(c)) as they contribute to the magnetic switching m = 5/2 —)■ 
m = —5/2. On the other hand the transitions m —)■ m + 1/2 compete for the total spin 
reversal and can be regarded as ‘backward’ processes (they follow the dashed lines in 
Fig. 2(c)). 

Furthermore, we distinguish two types of ‘forward’ transitions: (i) ‘absorption’ of 
spin-down electrons from the leads, i.e. the charging of the molecular orbital along the 
transitions |0,m) —?■ |l,m — 1/2)^ (full red arrow in Fig. 2(c)) and (ii) tunneling of 
spin-up electrons from the molecular orbital, i.e. a depletion process associated to the 
transitions |l,m — 1/2)^ —)■ |0,m — 1) (full blue arrow). Similarly one dehnes charging 
and discharging ‘backward’ processes (associated with the dashed lines in Fig. 2(c)). 
We hnd this analysis useful as it provides hints for a write-and-read scheme of states 
with well dehned molecular spin |0, m) when operating the SMM in the turnstile regime. 
Such a protocol will be discussed in the next subsection. 

3.2. The turnstile protocol 

A turnstile pumping cycle entails two steps: (i) the charging of the molecular orbital 
from the left lead while the drain contact is closed ixiif) 7^ 0) Xhif) = 0) and (ii) 
discharging/depletion through the drain lead (xiit) = 0, Xnit) 7^ 0)- We simulate the 
turnstile operation by appropriately tailoring the switching functions xl and xr fhe 
tunneling Hamiltonian Ht- 

The main idea behind the proposed operation is the following: use the charging 
sequences to prepare intermediate one-electron states via ‘forward’ tunneling from the 
left lead and write ‘empty’ molecular states |0,m) along ‘forward’ depletions to the 
right lead. To be more precise, let us discuss a single turnstile cycle at E = 0 in the 
NF conhguration, starting from the initial state |0,5/2). The associated transitions are 
depicted in Fig. 2(c). By opening the source (left) contact, the states |1,2)^ become 
populated by absorbing one spin - f electron (forward tunneling). In turn, if a spin 
- f electron is absorbed then the rightmost state |1,3) becomes activated (backward 
tunneling). Then, the left contact closes and the drain (right) electrode comes into 
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play. Now, the orbital is depleted through forward tunneling |1,3) |0,5/2) and 

|1,2)^ |0, 3/2), as spin - electron tunnels out into the right lead. An accurate 

operation would lead to the preparation of a single EMS or to an equally weighted 
mixture of EMS, but this scenario is not expected to work if the transverse anisotropy 
induces strong mixing of states |0,m). 

In view of this analysis, let us now discuss how this picture gets modihed in the 
presence of the transverse anisotropy. We start by describing the construction of the 
empty molecular states {| 930 , 1 ^) }• We hnd that if gfisB D the mixing of ‘empty’ 
molecular states |0,m) is negligible since the QTM between the states |0, =f 1/2) and 
|0, ±3/2) is very weak (see the dotted lines in Fig. 2(b)). In this case we hnd a one to 
one correspondence between v m as for any v in Eq.(lO) one can hnd a single m such 
that \^o,v) ^ |0,m). 

This simple correspondence fails as the magnetic held approaches resonant value 
gp^sBres = —D{Sz ± A'), and the Landau-Zenner tunneling processes between £^o,m and 
£o^rn'=m ±2 become important and lead to strong mixing of the states. Such a resonant 
regime will not be discussed in the present work. 

We now turn to Q = 1 states. For i? = 0, the states |l,m)=*= and |1,—m)=*= in 
Eq.(ll) are mixed by the transverse anisotropy term as and Sf_m are degenerate 
(see Eq.(8)). The mixing is indicated by the double arrow in Fig. 2(b). However, even 
a small magnetic held lifts this degeneracy and one hnds a rather small mixing of the 
states |l,m)^ and |1,—m)^ for E ^ 0. Once again, for each there is a single 

state |l,m')^ of whose weight |c®^,p in Eq.(ll) is by far the largest one. Under 
these conditions the one-to-one correspondence between the index n and the quantum 
number m is preserved for all states and allows us to index them as: 

Wpu) \^o,u) ^\Q,rn). (18) 

We shall close this section by noticing that although this representation works, v will 
always be read as an index, and not as a quantum number. Consequently, in the Iv^i^) 
basis the populations of the states will be denoted by and P\ipo^^)- 

3.3. Writing and reading the excited molecular states 

We performed transport calculations starting from the initial state |</?o, 5 / 2)5 so the density 
matrix describing the system at t = 0 is p(t = 0) = |<yCo, 5 / 2 )( 9 ?o, 5 / 2 |- As stated previously, 
we shall present results for small values of the ratio E/D ~ 10“^. The evolution of 
the relevant populations along a single turnstile cycle in the normal-ferromagnetic (NF) 
conhguration for E/D = 1/250 is presented in Fig. 3(a) and (b). The tunneling processes 
are similar to the ones discussed along Fig. 2(c) when F = 0. The state |<y 2 i, 3 ) is half 
hlled through spin-up ‘backwards’ tunneling whereas P\^p+^) + P\^p-.p} = 1/2- The small 
imbalance population of the states is due to the hnite J, while P\^+^) = P\ip-.p) 

at J = 0. Along the charging transition towards the Q = 1 sector, the population 
^10, 5 / 2 ) drops quickly to zero. The depletion cycle t G [2,4] ns simultaneously activates 
the states |<yCo, 5 / 2 ) and | 9 ?o, 3 / 2 ), the stationary regime being described by the RDO 
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Figure 3. (Color online) The evolution of the relevant populations along the first 
turnstile cycle for NF an AP configurations. The charging sequence corresponds to 
t € [0,2] ns and the depletion sequence to t G [2,4] ns. (a) and (b) - normal- 
ferromagnetic (NF) configuration, (c) - antiparallel (AP) configuration; the system 
ends up in a single excited state |</5o,3/2) (see also the discussion in the text). Other 
parameters: e = 0.25meV, J = O.lmeV, U = ImeV, ul = ImeV, ga = —ImeV, 
D = O.OdmeV, ggsB =0.005meV and E/D = 1/250. 


P = i\Pop/ 2 ){po, 5 / 2 \ + |¥^0,3/2)(<yi^0,3/2|)/2. Therefore one can nse the NF conhgnration 
to prepare an eqnally weighted mixtnre of states. Along the hrst depletion seqnence, in 
the NF conhgnration, (S'*) = 2 (see Fig. 4(a)). 

In the AP conhgnration, the hrst tnrnstile cycle drives the SMM ont of the gronnd 
state \po,b/ 2 ) directly into the hrst excited ‘empty’ molecnlar state |</?o, 3 / 2 ) (see Fig. 3(c)), 
with no fnrther mixing as in NF conhgnration. The popnlation of this excited state 
attains its maximnm valne within the depletion cycle. In the AP conhgnration the two 
EMSs can be viewed as binary digits (i.e. |v?o, 5 / 2 ) |0) and |(^o, 3 / 2 ) |1)) that are 

switched along the tnrnstile protocol. 

The writing of a single EMS depends crucially on the spin polarization of the leads. 
Fig. 3(c) indicates that the ‘backward’ processes (both charging and relaxation) are 
forbidden, as the state |(^i, 3 ) is not available along the charging sequence because the 
left lead provides no spin-up electrons. For the same reason there are no transitions 
from \ipi 2 ) back to |</?o, 5 / 2 ) on the discharging sequence. In fact, on the charging cycle 
the system occupies only two states with imbalanced populations P\^+^) > P\^-^) 
which eventually deplete in favor of |(po, 3 / 2 )- It is important to observe that on the 
depletion cycle the average total spin (S'*) = 3/2 and coincides with the molecular spin 
of the |(^o, 3 / 2 ) EMS (see Fig. 4(a)). 



























Quantum turnstile operation of single-molecule magnets 


13 


The reverse magnetic switching can be also implemented by simply reversing the 
bias {fiL ^ hi?) while keeping both contacts closed and then repeating the turnstile 
operation with the new initial state |<ho, 3 / 2 )- Then the system returns to |(po, 5 / 2 )- This is 
a classicalNOT operation, as the system evolves from |(po, 5 / 2 ) to |<ho, 3 / 2 ) and then back 
to |<ho, 5 / 2 ) without passing through a superposition of these states. 

Let us emphasize that the preparation of a pure excited molecular state |<ho,m) 
cannot be achieved in the standard transport regime. In that case the charge flows 
simultaneously to and from the leads and one cannot completely deplete the molecule 
and therefore {^o,m\pit)Wo,m) < 1- If should be mentioned that for larger S the time 
needed to achieve the full magnetic switching m = S ^ —S also increases as the system 
must visit all the intermediate states [26]. This fact suggests that the pair of consecutive 
states (|(po, 5 ), l^^o.s-i)) might be more appropriate for faster manipulation of magnetic 
qubits. 

Given these results one can ask about the time evolution of the total spin under 
repeated pumping cycles and on the possibility to read the states prepared along the 
turnstile operation by measuring currents. Fig. 4 summarizes our main results on 
transient currents and spin evolution along few turnstile cycles for the AP and NF 
conhgurations. The time-dependent occupation of the molecular orbital (the blue line 
in Fig. 4(a)) has a typical charging/relaxation pattern, with quick orbital hlling and 
slightly slower depletion. This can be seen by comparing the abrupt increase (in less 
than 1/2 ns) of the population at the beginning of the charging cycles (e.g. t = 4,8, 12ns) 
to the smooth tail of discharging which extends over Ins (e.g the time range [6, 7]ns). 

The total spin average (S'*) presented in Fig. 4(a) displays a step-like structure in 
both conhgurations. The steps scan both integer and half-integer values of (S'*), the last 
step for the AP conhguration corresponding to (S'*) = 3 being reached after t ~ 18 ns 
(not shown). In the AP conhguration the onset of half-integer steps corresponds to the 
depletion of the molecular orbital (Q = 1 —Q = 0), whereas the transition between 
half-integer to integer steps is associated to the charging process (Q = 0 —?■ Q = 1). 

The NF conhguration presents diherent features: there are fewer but longer steps 
of (S*), but no clear correspondence can be made between these steps and the behavior 
of the total charge Q. One notices that in this case the integer steps extend on some 
charging sequences and that half-integer values are encountered even if the orbital is 
empty. Fig. 4(b) displays the expected series of spikes for the transient currents in 
the AP conhguration. The period of the pumping cycles must be chosen appropriately 
in order to ensure full charging and discharging of the molecular orbital (we hnd the 
minimal period to be ~ 1 ns). 

We have found a similar behavior (not shown here) for the transient currents in 
the NF conhguration. The input current Jl vanishes when the orbital is fully occupied 
((5 = 1), whereas on the discharging sequence Jr drops to zero as the orbital depletes. 
The amplitudes of Jr and Jr are diherent because the charging process is faster than 
the depletion (see Fig. 4(a)). By inspecting Figs. 4(a) and (b) one observes that in the 
AP conhguration we have a one to one correspondence between the average (S'*) and the 
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Time (ns) 


Figure 4. (Color online) (a) The average total molecular spin (5*) in the anti¬ 
parallel (AP) and normal-ferromagnetic (NF) configurations - black lines. The step-like 
structure is discussed in the text. The total charge Q accumulated on the molecular 
orbital - blue line, (b) The transient currents Jl.r in the AP configuration. The 
half-integer steps of (Si) suggest that in the corresponding time range the SMM 
state is simply |0,m) - see the discussion in the text. The pumping period is 2ns. 
Other parameters: ul = 1 meV, fiR = —ImeV, e = 0.25 meV, J = 0.1 meV, 
U = ImeV and r = 0.5 meV, D = 0.04 meV, E/D = 1/250, ggsB = 0.005 meV, 
= V^ = 0.045 meV, hgT = 0.001 meV. 


peak-to-peak sequence in the transient currents: (S'*) acquires half-integer values only 
between a depletion peak and the next charging peak (e.g. for t G [2,4]ns (S'*) = 3/2), 
while between a charging peak and the next depletion peak the average spin is an 
integer. This means that the AP conhguration can be used to record experimentally 
the initialization of a given ‘empty’ molecular state \^o,m)- 

To this end it is sufficient to know the initial state of the molecule and to carefully 
‘count’ the transient peaks of Jl and Jr. We need to keep in mind though that 
Fig. 4(a) shows the average value of the total spin, which does not guarantee that 
along half-integer steps of (S'*) the system is in a pure state characterized by the RDO 
p ~ |0,m)(0,m|, especially for larger values of the transverse anisotropy when one 
expects stronger mixing of states. 
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Figure 5. (Color online) The populations P\ipo m) empty molecular states 

in the antiparallel configuration. Black line - 3 / 2)1 red line - ^^ 2)1 green line 

- ^lvo,-i/ 2 >i blue line - F|^ 3 __ 3 / 2 >i magenta line - (a) E/D = 1/250. (b) 

E/D = 1/75. (c) E/D = 1/25. On each turnstile cycle we indicate the dominant 
EMS. (d) P |/,3 _ 3 ) for different values of the ratio E/D. The other parameters are the 
same with those in Fig. 3. 

To further investigate the role of the anisotropy, we calculated the populations 
^Wo.Tn) of several ‘empty’ molecular states |v5o,m) for different values of the ratio E/D, 
at hxed magnetic held. Fig. 5(a) conhrms that at E/D = 1/250 the fc-th depletion cycle 
is well described by a single state \(po,m=b/ 2 -k)■ This proves the stepwise all-electrical 
writing of EMS (i.e. point (i) in the Introduction). 

By increasing the transverse anisotropy constant such that E/D = 1/75 we notice 
in Fig. 5(b) the emergence of a 2nd EMS on the depletion cycles. Nevertheless, the 
population of the dominant ‘empty’ molecular state exceeds 0.9 so we can still associate 
a well dehned molecular state to each of the depletion cycle. This no longer holds for 
E/D = 1/25. Fig. 5(c) reveals that the weight of the state |</?o,- 3 / 2 ) and |</?o,- 5 / 2 ) on 
the 2nd and 3rd depletion cycle increases up to 0.25, reducing the efficiency of the 
quantum turnstile protocol. Moreover, one can easily see that along the 4th cycle 
(i.e t G [14,16]ns) P\ipo_;i/ 2 ) + P\‘fioi/ 2 ) < f which suggests that other states have to 
be populated. We have found that the state |(pi^_ 3 ), corresponding to Q = 1 gets 
populated after the third cycle due to the forward transition |<po,- 5 / 2 ) —t |<yCi,- 3 ) via spin 
- I tunneling into the SMM. 
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Figure 6. (Color online) (a) The populations of active states along the second 
turnstile cycle for two values of the anisotropy constant. The dashed lines correspond 
to E/D = 1/75 and the solid lines to an increased value E/D = 1/25. More discussion 
is given in the text. Other parameters are as in Fig. 3. (b) and (c): Schematic 

representation of the tunneling processes between the states jt/sf and EMSs along 
the charging and discharging sequences. The numbers denote the dominant total 
magnetic moment m of the fully interacting one-particle states. 


In Fig. 5(d) we show that the population of this state becomes relevant with 
increasing the anisotropy, i.e. from a population of 0.03 at E/D = 1/100 to 0.25 
at E/D = 1/25. 

In order to explain the coexistence of two EMSs on the same depletion sequence 
when the transverse anisotropy increases we have to analyze the QTM between 
nearly degenerate Q = 1 states. By looking at the off-diagonal matrix element 
^(1, IKS*^ -b <S'^)|1, —1)^ we infer that by increasing E the hybridization of |1,1)^ and 
|1, —1)^ in the fully interacting states increases as well. We hnd that the weight of the 
‘minority’ state |1,—1)^ in increases from ~ 10“^ for E/D = 1/250 to ~ 10“^ 

for E/D = 1/25. The mixing of the states |1,2)^ and |1,—2)^ arises to the second 
order in E and is still negligible. As a consequence the accuracy of the hrst turnstile 
cycle is preserved even for E/D as large as E/D ~ 1/25 and that ~ This is 

conhrmed by the results presented in Fig. 5(a)-(c). In order to recover ‘clean’ EMSs on 
each depletion cycle for larger values of E/D one could slightly increase the magnetic 
held. The latter lifts even more the degeneracy of the states |1,1)^ and |1, — 1)^ and 
reduces therefore their mixing in the presence of E/D. 

Fig. 6(a) presents the relevant populations of EMSs and Q = 1 states on the 
second turnstile cycle {t G [4, 8]ns) for two values of the ratio E/D. For simplicity we 
plot the total population of states corresponding to the same dominant value of spin 
m. Figs. 6(b) and (c) indicate schematically the relevant tunneling processes between 
states of Hm along the charging and discharging sequences. The states are 

simultaneously hlled along the charging sequence. On the other hand, the states 
are less responsive to charging (this processes correspond to the dashed line in Fig. 6(b)) 
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because spin down tunneling is allowed only through the states whose weights 

are small. By similar arguments, one can see that the discharging process activates two 
EMSs, namely |(po,i/ 2 ) and |(po- 3 / 2 )- 

Both of these states acquire important weights in the RDO for E/D ~ 1/25 so the 
average total spin can no longer be associated to a well dehned value of the molecular 
spin. We therefore conclude that the enhanced QTM between Q = 1 states damages 
the efficiency of the turnstile protocol even if the EMSs involved in transport are not 
mixed. 

Let us note that the possibility to prepare a single EMS is not obvious as an open 
system is generally described by a mixed state. Our simulations also show that in the 
quantum turnstile regime one controls the transitions between any pair of intermediate 
molecular states (|(po,m), |<^o,m-i)) along a pumping cycle, in contrast to the full magnetic 
switching which involves only the pair (|v5o,s), |‘^o,-s))- 

Finally, we mention that if the SMM has an integer spin one cannot expect an 
accurate turnstile operation because the transverse anisotropy induces strong mixing 
between quasidegenerate EMSs (e.g. between |(po,i and |<po,-i)- 

4. Conclusions 

In the present work we address the transient transport regime and turnstile pumping 
across a single-molecule magnet coupled to external leads. The time-dependent 
evolution of the molecular states has been discussed in detail and signatures of the 
electrically induced magnetic switching on the transient currents were predicted. For 
ferromagnetic leads with antiparallel spin polarizations the turnstile protocol allows the 
stepwise writing and reading of excited molecular states. 

The evolution of the states along the turnstile operation can be ‘read’ indirectly 
from the behavior of the transient currents. More precisely, by recording the charging 
and discharging currents one can monitor the evolution of the system and identify the 
regimes where its density matrix is described by a single empty molecular state. This 
is somehow in contrast to the situation when the leads are simply normal metals and 
where the control of the excited spin states cannot be achieved as the molecular spin is 
reversed continuously. 

We show that the transverse anisotropy leads to the hybridization of nearly 
degenerate one-particle states which subsequently relax to empty molecular states with 
different values of the total spin. However, this dephasing effect can be reduced by 
applying a moderate perpendicular magnetic held. 

Another useful application of the turnstile regime that we address here is the 
possibility to mix several excited spin states (viewed as magnetic qubits) during the 
discharging cycles when the source electrode is normal and the drain electrode is 
ferromagnetic. Note that the short rise time of the switching functions used in our 
simulations is not essential for the turnstile operation and slower switching functions 
could be in principle selected to achieve a better resolution of the transient peaks. Our 
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predictive simulations clearly emphasize the potential of the molecular turnstiles as 
promising candidates for molecular spintronics. As a method approach we have used 
the generalized master equation formalism adapted to the turnstile configuration. 
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